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SUMMARY 


In  conducting  multiwavelength  scintillation  experiments  over  a  one 
mile  path,  it  was  found  that  a  finite-sized  transmitter  aperture  does  not 
represent  a  true  spherical -wave  source,  regardless  of  the  beam  divergence 
used.  This  observation  has  important  implications  for  the  interpretation 
of  other  data  on  the  "saturation"  of  scintillations,  and  in  the  present  experi¬ 
ments  has  lead  to  the  redesign  of  the  transmitter  to  achieve  a  virtual  point 
source.  This  has  been  completed,  including  a  provision  for  independently 
increasing  the  aperture  at  each  wavelength  up  to  two  inches  in  diameter,  for 
later  transmitter-aperture-dependence  measurements. 

A  highly-developed  thermal  probe  system  for  the  detailed  study  of 
turbulence  structure  was  made  completely  operational  f  and  a  portable  laser 
system  was  fabricated  for  independent  "optical  turbulence  strength"  measure¬ 
ments.  All  subsystems  of  the  experimental  facility  are  operating  without 
problems,  and  no  further  modifications  are  contemplated. 

In  order  to  put  the  problem  of  atmospheric  scintillations  in  current 
perspective  and  to  guide  the  comprehensive  experiments,  a  discussion  of 
recent  experimental  and  theoretical  issues  is  given.  It  is  felt  that  the  under¬ 
standing  of  the  turbulence  scattering  problem  has  progressed  to  the  point 
where  the  correct  questions  may  be  posed  and  definitively  investigated 
through  these  experiments. 
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I.  INTRODUCTION 

The  fundamental  goals  and  experimental  approach  for  this  program 
were  reviewed  in  the  preceding  Quarterly  Progress  Report.  During  the 
subsequent  period,  the  influence  of  the  transmitter  aperture  on  experimental 
results  was  found  to  be  substantial,  a  result  which  has  important  implica¬ 
tions  regarding  the  interpretation  of  "saturation"  of  scintillations  as  observed 
by  other  investigators.  The  multiwavelength  transmitter  was  modified  in 
order  to  achieve  a  definite  equivalent  to  a  point  source,  and  consequently, 
true  spherical -wave  saturation  measurements  may  now  be  conducted  for  the 
first  time.  In  addition,  separately  and  continuously  variable  apertures  are 
available  up  to  two  inches  in  diameter  at  each  transmitter  wavelength. 

A  highly -developed  thermal  probe  system  for  the  detailed  study  of 
turbulence  structure  was  made  completely  operational,  and  a  portable  laser 
system  was  fabricated  for  short-range,  "optical  turbulence  strength"  measure¬ 
ments.  All  subsystems  of  the  experimental  facility  are  operating  without 
problems,  and  no  further  modifications  are  contemplated. 

Highly  inclement  weather  conditions ,  including  consistent  fog  or 
rain,  precluded  the  gathering  of  comprehensive  data  subsequent  to  the  comple¬ 
tion  of  final  system  additions  described  above.  The  remainder  of  the  program 
will  be  devoted  entirely  to  the  use  of  the  system  to  obtain  data  relevant  to 
the  atmospheric  model,  the  wavelength  and  parameter -dependence  of  the 
onset  of  multiple  scattering  or  saturation,  and  the  influence  of  a  finite  trans¬ 
mitter  aperture  under  saturated  and  non-saturated  conditions.  As  discussed 
in  Section  III,  it  is  felt  that  the  understanding  of  the  turbulence  scattering 
problem  has  progressed  to  the  point  where  the  correct  questions  may  be 
posed  and  definitively  investigated  through  these  experiments. 

The  final  subsystem  developments  and  typical  experimental  results 
will  be  described  in  Section  II.  In  Section  III,  various  aspects  of  other 
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experimental  and  theoretical  efforts  will  be  discussed,  in  order  to  put  the 
overall  problem  in  current  perspective.  Future  plans  on  this  program  will 
be  restated  in  Section  IV. 

II.  PROGRESS  DURING  THIS  QUARTER 

A.  TRANSMITTER  APERTURE  DEPENDENCE 

Contrary  to  a  statement  in  the  preceding  report,  it  was  found 
experimentally  that  the  effects  of  a  finite,  diverging  transmitter  aperture 
are  highly  important,  and  that  the  realization  of  a  virtual  point  source  requires 
careful  design.  The  use  of  a  point  source  is  highly  desirable  when  the  funda¬ 
mental  aspects  of  saturat  on  are  being  explored,  i.e.  free  of  aperture  in¬ 
fluences,  and  the  implications  are  important  in  interpreting  the  saturation 
data  reported  by  other  investigators.  This  will  be  further  discussed  in 
Section  III. 


The  fact  that  divergent  but  nonzero  transmitter  apertures  fail 
to  yield  point-source  results  was  suggested  to  the  present  investigator  by 
Carlson^.  The  influence  of  the  transmitter  aperture  and  focusing  on  scintil¬ 
lation  results  has  been  rigorously  analyzed  within  the  single-scattering  realm 

2,3,4  5 

by  Ishimaru  and  Carlson,  ’  ’  and  semi-quantitatively  by  Fried.  Un¬ 

fortunately,  the  work  of  Ishimaru  and  Carlson  has  not  been  reduced  to 

numerical  curves  over  a  significant  range  of  parameter  values.  In  particular, 

4 

for  a  collimated  beam  it  has  been  shown  that  as  the  aperture  size  is  reduced 
from  a  very  large  Fresnel  number  (N)  to  a  very  small  N,  the  variance  of 
log-amplitude  scintillations  should  decrease  from  the  plane-wave  value  through 
a  minimum  (at  N  =  1)  and  then  increase  to  the  final  value  for  a  spherical 
wave  or  point  source^  (Figure  1).  For  a  focused  beam,  the  minimum  is  even 
sharper,^  representing  "transmitter  aperture  averaging."  For  the  trans¬ 
mitter  as  described  in  the  preceding  Quarterly  Report,  the  value  of  N  at  the 

o 

shortest  (4880  A)  wavelength  was  on  the  order  of  unity,  so  that  a  change  to  a 
point  source  would  be  expected  to  result  in  an  increase  in  scintillation. 
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However,  as  shown  in  Figure  1,  the  opposite  effect  occurred 
in  our  experiments.  Starting  with  a  one -inch  diameter  beam  which  was 
diverged  to  a  1  milliradian  angle,  successively-smaller  aperture  stops  were 
placed  over  the  output  window.  The  results  are  summarized  in  Table  I,  and 
indicate  a  significant  decrease  in  scintillations  for  a  point  source  as  opposed 
to  a  source  with  N  ~  1.  The  change  was  less  pronounced  under  high-scintilla¬ 
tion  conditions  ,  but  in  either  case  is  readily  repeated  and  easily  discernible 
by  eye . 


Aperture  Size  Log-Amplitude  Variance 


1" 

0.  13 

0.45 

1/4" 

0.05 

0.32 

1/8" 

0.04 

0.25 

Table  I.  Variance 

vs.  Transmitter  Aperture 

(4880  A) 

(One  Mile  Path) 

The  theoretical  curve  in  Figure  1  fails  to  predict  the  experi¬ 
mental  results  because  the  transmitter  was  highly  divergent.  The  param¬ 
eter  a  2  ,  which  indicates  defocusing,  was  much  larger  than  explicitly  covered 
in  the  analysis  of  Ref.  3,  so  that  actual  numerical  predictions  for  this  situa¬ 
tion  have  not  yet  been  given. 

The  effect  observed  here  has  important  implications  in  rela¬ 
tion  to  other  commonly -quoted  measurements  of  saturation,  as  discussed 
further  in  Section  III.  For  the  present  case,  it  became  apparent  that  a 
diverging  beam  is  insufficient  for  the  application  of  point-source  theory. 

In  order  to  realize  a  multiwavelength  transmitter  with  independ¬ 
ent  aperture -adjustment  capabilities  at  each  wavelength,  including  the 
achievement  of  a  virtual  point  source,  it  was  necessary  to  eliminate  the  con¬ 
straint  represented  by  the  single  9  kHz  chopping  wheel  located  at  the  input 
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focus  of  the  transmitting  telescope.  As  reported  previously,  we  have  been 
able  to  eliminate  synchronous  detection  while  retaining  our  dynamic  range, 
and  hence  we  are  permitted  to  use  separate,  smaller  chopping  units  at  the 
foci  of  beam-forming  lenses  at  the  output  of  each  of  the  three  lasers.  Through 
the  use  of  electronic  etching  techniques  and  ultra -precise  machining,  choppers 
were  developed  which  are  free  of  residual  line-locked  vibration  and  wobble, 

hence  eliminating  small  instrumental  spectrum  and  covariance  components  in 

7 

the  received  scintillations. 

The  point  source  configuration  of  the  transmitter  was  designed 
by  first  choosing  a  diffraction-limited,  collimated  divergence  of  two  milli- 
radians  at  4880  A  and  10.6  microns,  which  represents  sufficient  spread  for  a 

g 

"spherical  wave"  while  preserving  a  high  enough  power  density  on  the  small 
receiver  apertures  to  maintain  the  required  dynamic  range.  Refracting  optics 
were  then  utilized  to  match  the  two  respective  lasers,  through  the  combining 

mirrors,  to  the  off-axis  parabolic  output  telescope  such  that  beam  waists  of 

9 

the  proper  size  appear  at  the  output  window.  Furthermore,  by  sliding 
certain  lenses  beyond  the  chopping  wheel,  the  output  beam  sizes  at  these  two 
wavelengths  may  be  independently  varied  up  to  two  inches  in  diameter  (un¬ 
obstructed).  In  the  point-source  condition,  the  inverse  Fresnel  number 
(a  i  L  in  Refs .  2-4)  is  4 , 200  at  4880  A  and  194  at  10.6  microns  ,  resulting  in 
definite  point-source  characteristics. 

At  the  intermediate  wavelength  of  1.15  microns,  the  laser 
power,  mode  limitations,  and  dynamic  range  conditions  are  such  that  the 
collimated,  two-inch-diameter  beam  cannot  be  decreased  in  size,  since  greater 
divergence  would  result.  However,  since  the  laser  radiation  has  substantial 
transverse  mode  structure,  the  effective  aperture  may  be  expected  to  be 
significantly  smaller  than  the  physical  beam  size. 


^Further  aperturing  at  10.6  microns  showed  no  effect  on  measured  variances. 
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B. 


MEASUREMENTS  WITH  POINT-SOURCE  TRANSMITTER 


An  example  of  simultaneous ,  multiwavelength  data  gathered 
with  the  point-source  configuration  at  4880  A  and  10.6  microns  is  given  in 
Table  II.  The  data  were  taken  in  the  earl/  afternoon,  with  hazy  sky  and  light 
wind  conditions.  The  probability  distributions  at  each  wavelength  were 
accurately  log  normal. 


4880  A 

1.  15p 

10. 6  p 

Log -amplitude  variance  from 
computer 

0.36 

0.26 

0.011 

Log -amplitude  variance  from 
probability  slope 

0.38 

0. 24 

0.012 

Amplitude  correlation  length 
at  receiver  (cm) 

2.38 

- 

11.0 

Bandwidth  of  scintillations  (Hz) 

35 

22 

Table  II.  Scintillation  Data  with  New  Transmitter 
Configuration  (One  Mile  Path) 


The  ratio  of  variance  at  10.6  microns  vs.  4880  A  is  33,  which 
may  be  compared  to  the  theoretical  (k  7^6 )  value  of  36.  However,  the  variance 
at  1.15  microns  is  70%  higher  than  its  theoretical  value,  relative  to  4880  A, 
thus  indicating  that  transmitter  aperture  effects  are  important  at  that  wave¬ 
length;  this  is  unavoidable,  as  discussed  above. 

The  amplitude  correlation  length  of  2.38  cm  at  4880  A  may  be 

1/ 

compared  with  the  theoretical  (X  L)/  value  of  2.80  cm,  while  the  experi¬ 
mental  and  theoretical  values  at  10.6  microns  are  11.0  and  12.9  respectively. 
Due  to  a  temporary  malfunction,  no  value  was  obtained  at  1.15  microns. 
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The  scintillation  bandwidth*  were  quite  small,  due  to  the  light 
wind  conditions.  The  ratio  of  bandwidth*  at  4880  A  and  1.15  microna  was 
1.59,  compared  to  a  theoretical  (\)  ratio  of  1.54,  thus  indicating  that  the 
covariance  may  not  be  seriously  affected  by  the  finite  transmitter  at  1.15 
microns.  The  bandwidth  at  10.6  microns  was  too  low  for  measurement  under 
such  light  wind  conditions;  thiu  is  not  the  usual  case. 

It  is  not  clear  from  the  above  data  whether  the  variances  were 
saturated  during  this  run.  This  is  discussed  further  in  Section  III.  As  soon 
as  weather  permits,  comprehensive  data  will  be  taken  to  establish  the  true 
point. source  wavelength-dependence  of  saturated  variance  and  covariance, 
as  discussed  in  the  previous  report  and  Section  III.  In  particular,  the  previously 
mentioned  possibility  that  saturated  variances  follow  the  k  law  will  be  proven 
or  disproven. 


It  should  be  noted  that  under  conditions  of  substantial  rain, 
which  has  prevailed  recently,  the  scintillations  are  observed  to  be  much 
faster,  with  spectra  at  each  wavelength  up  to  500  Hz,  and  very  small  ampli¬ 
tude  correlation  lengths.  The  effects  of  turbulence  scattering  cannot  be 
studied  under  such  conditions. 

C.  OPTICAL  TURBULENCE  STRENGTH  MEASUREMENT 
SYSTEM 

It  was  decided  that  there  would  be  substantial  value  in  incor¬ 
porating  a  provision  for  measuring  the  "optical  strength  of  turbulence"  using 
a  short-range  laser  link,  as  described  in  Ref.  10.  Since  the  scintillation 
from  a  laser  at  e.g.  500  feet  is  small,  the  single-scattering  (Rytov)  theory 
applies  without  questions  of  saturation,  and  an  independent  determination  of 
the  refractive  index  structure  constant  (C  )  is  obtained.  ^ 
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In  order  to  achieve  this,  we  have  modified  1-milliwatt  He-Ne 
laser  (6328  A)  to  incorporate  electronic  chopping,  and  have  mounted  the  unit 
on  a  tripod  and  cart  with  portable  filament  and  high-voltage  better/  supplies. 
The  laser  is  bore-sighted  with  a  rifle  scope,  and  may  be  steered-in  rapidly. 
The  beam  is  shape  to  comprise  a  true  point  source.  The  receiver  is  equipped 
with  an  S-20  photomultiplier  and  10-angstrom  background  filter  for  this  signal, 
and  the  variance  is  rapidly  measured  using  the  regular  system  computer. 

This  system  is  complementary  to  the  thermal  instrument  which  will  now 
be  described. 

D.  THERMAL  MICROPROBE 

During  the  past  quarter,  a  highly  sensitive  and  stable  thermal 

microprobe  for  the  direct  measurement  of  the  refractive  index  (temperature) 

structure  constant  became  operational.  The  instrument  is  a  refinement 

12-14 

of  a  concept  previously  described  ,  and  incorporates  feedback-stabiliza¬ 

tion  against  drift  and  a  low-noise  subcarrier  bridge  technique.  The  elec¬ 
tronics  and  operational  characteristics  are  discussed  in  detail  in  Appendix  A, 
along  with  calibration  and  data -reduction  procedures.  The  capabilities  of  the 
instrumentation  will  be  reviewed  in  the  present  section. 

The  instrument,  shown  in  Figure  2,  consists  of  two  fine  wire 

thermal  probes  with  a  variable  separation;  and  portable,  tripod -mounted 

electronics  which  process  the  instantaneous,  fluctuating  difference  in  the 

temperatures  at  the  probes.  The  mean-square  of  this  difference  is  presented 

on  a  built-in  meter  with  a  long  averaging  time  constant.  Appropriate  gain 

settings  are  chosen  with  the  aid  of  a  fast-response  saturation  indicator.  The 

dynamic  range  of  the  instrument  covers  C  2  values  from  1. 2  x  10  11  m"  ^  3 

n 

to  6  x  10  17  m~2/3,  which  more  than  brackets  real  values  at  the  ground. 

An  output  which  is  proportional  to  the  instantaneous  tempera¬ 
ture  fluctuations  is  provided  which,  through  the  use  of  a  long  cable,  permits 
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the  spectral  and  probability  distribution  analysis  of  the  fluctuations.  These 
analyses  utilize  the  scintillation -computer  electronics. 

An  example  of  temperature  fluctuations  is  shown  in  Figure  3. 
The  vertical  scale  calibration  is  0.07°C  per  division  and  the  sweep  speed 
is  500  msec/cm.  On  this  particular  day,  the  thermal  probe  indicated  a  C^2 
of  2.9  x  10  15  m  2/3,  and  the  wind  was  very  light.  On  highly  turbulent  days, 
fluctuations  of  more  than  10  times  those  shown  in  Figure  3  may  be  observed. 

The  probability  distribution  of  fluctuations  such  as  those  shown 

in  Figure  3  is  expected  to  be  gaussian.  It  must  be  noted,  however,  that 

10 

recent  experiments  have  indicated  the  presence  of  anomalous  spikes  which 
may  indicate  plumes  or  other  irregularities  near  the  ground.  We  have  not 
seen  such  irregularities,  nor  the  attendant  "tail"  on  the  probability  distribu¬ 
tion,  but  weather  conditions  have  permitted  little  data -taking  since  the  system 
became  operational.  Under  stable  conditions,  the  probability  tail  may  be 
expanded  and  examined  in  detail. 

Under  these  conditions  of  light  wind,  it  is  not  unusual  for  the 
atmospheric  statistics  to  be  highly  nonstationary  over  time  periods  neces¬ 
sary  for  properly -averaged  measurements.  For  example,  immediately  sub¬ 
sequent  to  the  meter  readout  of  C^2  =  2.9  x  10  15  m  2/3  above,  probability- 
slope  determinations  of  <(Ai)2>  showed  values  which  changed  abruptly  and 
ranged  from  4.5  x  10  15  to  7.8  x  10  16.  An  example  is  shown  in  Figure  4,  in 
which  CR2  changed  suddenly  during  a  probability  run.  The  points  are  tran¬ 
scribed  to  probability  paper  from  the  continuous  ,  smoothed  experimental 
plot  obtained  on  an  XY  recorder. 

The  spectrum  of  the  fluctuations  of  Figure  3  is  shown  in 
Figure  5.  For  an  average  wind  velocity  of  0.5  m/sec,  the  observed  cutoff 
frequency  of  60  Hz  implies  an  inner  scale  of  turbulence**  of  0.75  cm.  This 
rather  large  value  indicates  the  low  energy  at  such  low  wind  levels,  and  can 
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appreciably  affect  the  scintillation  statistics,  as  discussed  further  in 
Section  III.  As  discussed  in  Appendix  A,  it  is  possible  to  investigate  the 
theoretical  f"  5^6  dependence  of  the  rms  AT  spectrum  using  the  present  equip¬ 
ment  and  a  slightly  modified  technique. 

Within  the  inertial  subrange,  the  Kolmogorov  theory**  predicts 
that  the  mean-square  temperature  difference  between  the  two  thermal  probes 
will  vary  as  the  2/3  power  of  their  separation.  This  will  be  investigated 
experimentally  with  the  above  system. 

A  preliminary  measurement  of  the  height -dependence  of  C  2 
near  the  ground  showed  the  value  to  be  40%  lower  at  81  inches  than  at  45  inches; 
hence,  the  determination  should  be  made  at  the  actual  beam  height.  In  addi¬ 
tion,  under  stable  conditions  the  uniformity  along  the  path  may  be  checked, 
utilizing  the  portability  of  the  instrument.  However,  due  to  the  coincident 
path  nature  of  these  multiwavelength  experiments  ,  this  uniformity  is  not  highly 
crucial. 


As  of  this  writing,  there  has  not  yet  been  an  opportunity  to 
compare  the  optical  and  thermal  with  the  three -wavelength  variances  of 
log  amplitude  scintillation.  However,  on  one  occasion  the  thermal  was 
measured  simultaneously  with  the  one-mile  scintillations  at  4880  A  on  a  low- 
turbulence  day,  resulting  in  a  probably-fortuitous  agreement  of  7%  of  the 
theoretical  variance. 

To  summarize,  development  is  complete  on  these  instruments, 
which  will  now  be  utilized  as  follows: 

1.  All  simultaneous  multiwavelength  measurements  of 

log -amplitude  variance,  covariance,  probability  distribution ,  and  amplitude 

scintillation  will  be  backed  up  by  concurrent  thermal  and  short-range  optical 

determinations  of  C  2  . 

n 
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2.  Experiments  relative  to  the  atmospheric  model  per 
se  will  be  conducted  to 

a.  Investigate  any  non-gaussian  or  spiking 
behavior  in  the  temperature  fluctuations. 

b.  Investigate  the  validity  of  the  "2/3  law"  through¬ 
out  the  inertial  subrange  through  variable 
probe  -separations . 

c.  Investigate  the  extent  of  the  inertial  subrange, 
using  "b"  for  the  outer  scale  (and  possibly 
large  inner  scale  conditions)  and  spectral  cut¬ 
off  of  the  inner  scale. 

d.  Investigate  the  validity  of  the  theoretical  5/3 
(or  rms  5^6)  spectral  characteristic  of  thermal 
fluctuations . 

e.  Investigate  path  length  and  height  dependence 

of  C2 
n 

The  validity  and  extent  of  the  subrange  is  especially  of  interest  during  inverted 
and  low-wind  conditions,  as  discussed  in  Section  III. 

III.  DISCUSSION  OF  THE  CURRENT  STATE  OF  THE  TURBULENCE 
SCATTERING  PROBLEM 

As  a  result  of  discussions,  meeting,  and  preliminary  experimental 
results,  it  has  become  apparent  that  a  critical  review  of  recent  experiments 
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and  theoretical  efforts  is  desirable  in  order  to  put  the  overall  problem  in 
perspective  and  to  guide  the  experiments  which  are  now  commencing. 

In  discussing  those  aspects  of  scintillation  which  are  not  understood, 
we  must  address  ourselves  to  (1)  the  range  of  validity  of  the  atmospheric 
model,  (2)  the  mysterious  problem  of  the  parametric  dependence  of  satura¬ 
tion  or  multiple  scattering,  and  (3)  the  influence  of  transmitter  apertures. 
We  will  first  review  the  results  of  other  experimental  efforts,  and  then  con¬ 
sider  the  theoretical  status. 


A.  EXPERIMENTS 

15,16 

The  saturation  of  scintillations  has  been  observed  in  Russia, 

17  18  8  19  10 

and  by  Deitz,  ’  Lawrence  and  co-workers,  '  and  Fried  .  It  is 

important  to  review  just  what  was  measured,  and  how  well  the  results  present 

a  consistent  picture. 

As  discussed  in  Appendix  B,  there  are  several  second -order 
moments  which  have  been  used  by  experimentalists  and  theoreticians.  Some 
of  the  inter-relationships  (Appendix  B)  are  dependent  on  an  assumption  of 
log  normality  of  the  scintillations.  Furthermore,  ambiguities  and  actual 
errors  in  these  interrelationships  have  appeared  in  the  literature.  A  con¬ 
sistent  notation  and  definitions  are  given  in  the  Appendix  and  summarized  here: 


C/(0)  =  log  amplitude  variance 

C  (0)  =  log  intensity  variance  =  4Cj(0) 

ffj2  =  normalized  intensity  variance 

—  log^  <  I2>  /<I>2,  where  I  is  the  intensity 
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utilized 


In  the  early  saturation  measurements,  Gracheva^’  ^ 

an  incoherent  source  in  a  quasi-plane-wave  configuration  (80  mm  diameter 

lens  with  a  "slightly  divergent"  beam)  to  measure  <r,£  .  Then,  with  the 

assumption  of  log  normality,  the  saturation  of  Cj(0)  at  a  value  of  2.  56  is 

inferred.  DeWolf*^  has  pointed  out  an  error  in  the  data  reduction,  such  that 

the  saturated  Cj(0)  is  1/4  of  this  or  0.64,  implying  a  saturated  log  amplitude 

variance  Cf(0)  =0.16.  In  view  of  subsequent  measurements  by  others,  it 

must  be  concluded  that  this  value  is  too  low.  It  should  be  pointed  out  that 

the  source  was  not  well  defined,  and  that  the  assumption  of  log  normality  is 

21 

quite  critical  since  small  departures  can  result  in  an  apparent  saturation. 


Other  investigators  have  worked  with  quasi -point-sources 
(spherical  waves ) ,  as  follows: 

Deitz^  utilized  a  sophisticated  photographic  technique  to 
determine  o^2  and,  again  assuming  log  normality,  inferred  a  saturated  value 
of  0^(0)  which  essentially  agrees  with  that  of  the  corrected  Gracheva  data. 

In  this  measurement,  the  transmitter  had  an  approximately  1.7  mRad  diver¬ 
gence  and  5  cm  diameter,  and  it  was  found  that  a  divergence  decrease  to  one 
mRad  resulted  in  an  increase  in  the  saturated  variance  to  unity.  The  cor¬ 
responding  Cj.(0)  is  0.25.  In  view  of  the  discussion  of  Section  II,  it  is  apparent 
that  the  larger  divergence  was  not  sufficient  to  ensure  point-source  behavior, 
due  to  the  5  cm  aperture:  and,  in  fact,  the  general  conclusion  of  Figure  1, 
i.e.  that  divergent,  large  apertures  lead  to  greater  variances  than  point 

sources,  is  supported  in  Deitz's  non-saturated  data.  It  should  be  noted  that 

21 

the  influence  of  the  log-normality  assumption  was  again  a  factor. 

1 8 

In  more  recent  measurements  ,  Deitz  and  Alcaraz  utilized  a 
5  cm  transmitter  with  a  1.9  milliradian  divergence  to  obtain  saturation  at 
Cf(0)  =  0.25  with  a  fall-off  to  0.  16.  In  this  work,  the  log  normality  assump¬ 
tion  was  removed,  since  both  C^(0)  and  <rj  were  measured. 
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Lawrence  utilized  a  5  cm  transmitter  with  a  1  milliradian 
divergence  to  directly  obtain  a  saturated  C^(0)  of  0.7.  Again,  the  variances 
in  the  non-saturated  region  were  larger  than  predicted  by  theory,  as  in 
Figure  1 ,  and  the  transmitter  may  have  failed  to  represent  a  point  source. 

In  the  recent  Emerson  Dry  Lake  experiments,^  transmitter 

apertures  of  2-5  mm  were  employed.  At  3  mm,  the  inverse  Fresnal  number 

is  approximately  25,  so  that  these  measurements  may  have  represented  a 

better  approximation  to  a  point  source  determination.  C^(0)  was  again  directly 

measured  to  saturate  at  a  maximum  of  0.7.  The  theoretical  predictions  were 
5 

"corrected"  for  the  finite  aperture,  but  in  the  wrong  direction  (predicted 
variances  less  than  for  a  point  source).  This  is  again  a  manifestation  of  the 
phenomenon  of  Figure  1:  deliberate  divergence  results  in  greater  scintillation. 

It  is  seen  that  there  is  a  significant  range  in  the  "saturated" 
variances  reported  to  date.  The  higher  value  of  0.7  for  C^(0)  is  substantiated 
by  measurements  made  earlier  on  this  program  with  a  finite  aperture.  How¬ 
ever,  we  are  only  now  beginning  experiments  with  a  true  point  source  trans¬ 
mitter,  which  should  resolve  the  question.  The  sample  measurements  in 
Table  II  indicate  a  variance  at  visible  wavelengths  which  is  definitely  in  excess 
of  the  Deitz  and  Russian  values;  furthermore  since  independent  C*  determina¬ 
tions  were  not  included,  it  is  impossible  to  say  whether  or  not  this  variance 
yet  represents  saturation. 

The  importance  of  determining  the  variances  and  parameter 
breakpoint  for  saturation,  free  of  transmitter  aperture  influences,  is  now 
clear.  It  may  be  noted  that  log  amplitude  variances  of  0.25  and  0.7  represent 
greatly  different  scintillation-dynamic-ranges,  which  is  of  considerable 
practical  importance. 

B.  THEORETICAL  ASPECTS 

The  well  known  single-scattering  (Rytov)  prediction  for 
log-amplitude  variance  of  a  spherical  wave  is^ 
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(1) 


C,(0)  =  0.124  C2  k7/6  Lu/6 


n 


where  k  is  the  optical  wavenumber  and  L  is  the  distance  (range).  The  plane 

wave  prediction  is  2.5  times  larger.  ^  Also,  within  this  Rytov  realm,  solu- 

2-4 

tions  for  transmitter  aperture  effects  have  been  obtained.  Theoretical 
efforts  are  currently  concentrated  on  predicting  at  least  the  parameter  break¬ 
point  and  general  character  of  saturation  (multiple  scattering) , 
but  it  is  not  clear  that  definitive  progress  has  been  made,  as  discussed  below. 
The  important  question  of  the  consistency  and  range  of  validity  of  the  atmos¬ 
pheric  model  is  also  under  active  consideration,^  *  ^ »  1®  *  ^5  ancj  will  also  be 
discussed  in  this  section. 


1. 


Saturation 


Theoretical  predictions  of  saturation  involve  the  break¬ 
down  of  the  Rytov  approximation,  which  for  many  purposes  seems  no  more 
useful  than  the  first-Born  approximation.  There  are  inherent  inconsistencies 
in  the  Rytov  approximation  outside  the  first-Born  region.  In  addition,  there 
is  some  confusion  in  the  literature  over  the  nonzero  value  of  the  mean  log 
amplitude  and  the  distinction  between  variances  and  second  (noncentral) 
moments.  This  will  be  discussed  in  a  later  report. 

Theoretical  investigators  have  been  generally  con¬ 
cerned  with  expressing  an  actual  or  true  C.(0)  vs.  the  Rytov-predicted 

2  26 

value  of  Eq.  (1),  as  the  turbulence  (C  )  and  path  length  (L)  increase.  In 

20  " 

such  an  effort,  DeWolf  predicted  an  asymptotic  saturated  value  for  C^(0) 
of  In  2  as  the  Rytov-predicted  value  approaches  infinity.  This  is  consistent 

with  the  corrected  Russian  data.  However,  as  pointed  out  above,  it  is 

%«* 

v 

significantly  less  than  experimental  values.  The  analysis  furthermore 

*It  should  be  noted  that  experimental  results  for  C/tj(0)  vs.  ^  ^(0)  (in¬ 

creasing  Cn2  or  path  length)  show  a  maximum  and  small  negative  s^op'e  as  the 
Rytov-predicted  value  continues  to  increase,  so  that  the  asymptote  may  be 
smaller  than  the  values  observed  under  moderately  saturated  conditions. 
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predicted  Rayleigh,  rather  than  log  normal,  statistics  for  strong  scintilla¬ 
tions.  Like  the  Russian  experiments,  DeWolf's  analysis  was  in  terms  of 

or  <r2  ,  whose  relationship  to  C  .(0)  (as  measured  in  the  U.S.)  depends 

T  f  21 
rather  sensitively  on  the  log  normality  of  the  statistics. 

22 

In  more  recent  work,  DeWolf  has  acknowledged  the 
l3g  normal  --  rather  than  Rayleigh  --  nature  of  various  experimental  results. 

It  was  surmised  that  high-irradiance  fluctuations  are  sufficiently  separated  in 
time  that  each  involves  a  different  instantaneous  path,  which  leads  to  log  normal 
rather  than  Rayleigh  statistics.  He  then  assumes  that  the  departure  from 
Rayleigh  statistics  in  his  saturation  calculation  occurs  far  enough  out  on  the 
tail  of  the  distribution  to  leave  that  calculation  valid.  However,  in  a  final 
footnote  he  acknowledges  that  values  of  C^(0)  measured  in  the  U.S.  have  been 
higher  than  In  2,  so  that  even  this  assumption  is  questionable.  It  is  there¬ 
fore  apparent  that  much  more  work  is  required  here. 


In  another  recent  theoretical  effort,  the  renormaliza¬ 
tion  technique  and  ordering  arguments  have  been  applied  to  obtain  a  satura¬ 
tion  result  similar  to  but  somewhat  more  detailed  than  DeWolf's.  Although 
the  validity  is  unestablished,  the  analysis  has  very  interesting  aspects.  Un¬ 
fortunately,  this  analysis  to  date  has  employed  a  gaussian  rather  than 
Kolmogorov  atmospheric  model,  which  xeads  to  parameter  dependencies  which 
may  differ  seriously  from  those  of  the  "correct"  model.  For  instance,  the 
log  amplitude  variance  in  the  Rytov  realm  depends  on  (An)2  /  Lk2,  where  n 
is  the  refractive  index  and  t  is  a  turbulence  scale  factor.  This  is  sufficiently 
different  from  the  C^2  k7^6  L11^6  prediction  of  Eq.  (1)  to  render  any  resultant 
parametric  conclusions  doubtful.  The  analysis  is  again  in  terms  of  <r,£  , 
which  relates  to  the  measured  C.  T(0)  only  for  log  normal  statistics. 


infinity  is 


One  result  is  that  the  predicted  limit  as  (An7)  goes  to 
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lim  <r^. 


In  [  2  + 


] 


(2) 


4m2 
4m2  +  1 


where  m  =  — —  .  This  asymptote ,  which  may  be  identified  with  Cj(0)=  4  C^(0) 
for  log  normal  statistics,  ranges  between  In  3  and  In  2,  the  higher  value 
representing  short  wavelengths  and/or  path  lengths.  For  visible  wavelengths, 
the  predicted  value  is  still  low,  and  the  wavelength-dependence  (e.g.  at 
10  microns)  is  obviously  quite  weak,  which  may  or  may  not  agree  with  our 
future  experiments.  The  limit  for  constant  and  increasing  L  is  In  2,  as 
in  Refs.  20,22. 

It  is  interesting  to  note  that  the  complete  curves^  for 

cr2  vs.  C  (0)  exhibit  the  maximum  and  fall -off  reported  by 

T  predicted  I  Rytov 

experimenters.  The  general  results  involve  a  somewhat  complex  function, 
from  which  a  broad  range  of  numerical  curves  have  not  yet  been  derived. 
Agreement  with  the  Russian  data^5,  ^  is  demonstrated,  although  the  param¬ 
eter  l  is  chosen  to  effect  this  agreement.  Rayleigh  statistics  are  predicted 
for  the  limiting  case. 


It  is  to  be  hoped  that  the  approach  of  Ref.  23  will  be 
extended  to  the  Kolmogorov  atmosphere,  so  that  the  weak  wavelength -depen¬ 
dency  and  low  value  of  saturated  at  visible  wavelengths  may  be  critically 
examined.  If  our  experiments  firmly  establish  a  low  value  of  variance  satura¬ 
tion  at  10.6  microns,  the  existing  analysis  will  be  disproven;  on  the  other 

hand,  if  the  10. 6^  point-source  results  for  high  C2  continue  to  agree  with 

n 

Rytov  predictions  up  to  levels  comparable  to  saturation  in  the  visible,  the 
present  prediction  will  merit  considerable  attention. 

2 .  Transmitter  Aperture  Dependence 

Since  the  gross  details  of  saturation  with  a  point  source 
are  not  yet  understood,  it  is  clear  that  little  can  be  said  about  transmitter 
aperture  dependence  under  saturated  conditions  before  the  experiments  are 
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performed.  However,  as  pointed  out  in  Ref.  18,  some  of  the  details  of 
apparent  saturation  curves  may  involve  transmitter  aperture  effects  within 
the  single-scattering  realm.  Referring  to  the  theoretical  curve  in  Figure  1, 
if  we  start  out  in  the  near  field  of  a  collimated  beam  and  then  increase  the 
range  L,  the  value  of  N  increases  proportionately.  Hence,  the  log  amplitude 
variance  will  grow  more  slowly  than  L11^6,  and  a  resultant  curve  of  observed 
vs.  theoretical  point-source  (or  plane  wave)  variance  will  exhibit  some  of  the 
characteristics  of  saturation.  The  effect  is  not  strong  enough  to  explain  the 
reported  range -dependent  saturation  measurements,  but  can  definitely  affect 
the  knee  of  the  saturation  curves  in  Refs.  15-17. 

This  phenomenon  illustrates  the  importance  of  obtain¬ 
ing  fundamental  saturation  data  with  a  virtual  point  source.  At  the  same  time, 

experiments  to  confirm  the  transmitter  aperture  theory  under  single-scatter - 

1 8 

ing  conditions  are  highly  important. 


3 .  The  Atmospheric  Model 

The  range  of  validity  of  the  inertial  subrange  model 

has  yet  to  be  firmly  established.  On  very  calm  days,  for  instance,  there 

seems  to  be  insufficient  energy  to  establish  the  subrange,  and  hence  there 

is  no  applicable  theory  even  though  anomalously  large  and  non -log -normal 

10 

scintillations  are  observed.  Even  under  "normal"  conditions,  there  may 
be  a  non-gaussian  or  spiking  behavior  of  An  (Ref.  10),  which  would  negate 

24 

the  usual  argument  for  a  log  normal  (Rytov)  wave  function.  This  may  be 

27  28 

related  to  an  observed  breakdown  in  2/3 -law  behavior  near  the  ground.  ’ 

These  spikes  represent  refractive  index  spectral 
characteristics  outside  the  outer  scale,  and  necessitate  probability-slope 
deriviations  for  a  useful  .  Furthermore,  if  the  outer  scale  is  sufficiently 
small  that  substantial  portions  of  the  optical  filter  function  are  weighted  by 
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a  non -Kolmogorov  spectrum,  then  existing  theory  is  inadequate  for  even 
single-scattering  predictions.  This  condition  is  also  postulated  as  occurr¬ 
ing  under  diurnal  inversion  conditions,  during  which  qualitatively  different 

scintillation  patterns  have  been  observed  and  C  measurements  (low)  have 

n  17 

ceased  to  correlate  with  variance  measurements  (high).  A  further  con¬ 
firmation  of  this  condition  was  that  temperature  correlations  increased  with 

increasing  probe  separations,  indicating  that  the  separations  were  larger 

.  .  .  17 

than  the  outer  scale. 

Even  when  the  inertial  subrange  is  established  and 

the  outer  scale  is  "typical"  it  is  important  to  investigate  the  effects  of  a 

nonzero  inner  scale.  The  inner  scale  is  believed  to  be  typically  1  mm,  and 

is  assumed  to  be  zero  in  many  single-scattering  theoretical  developments. 

(The  usually-stated  condition  is  that  the  receiver  must  be  in  the  far  field  of 
6  1 1  4 

the  inner  scale.  *  )  Carlson  has  shown  that  an  inner  scale  on  the  order 

of  1  cm  or  more  can  substantially  reduce  scintillations  for  point-source 

transmitters.  The  analysis  simply  involves  setting  the  refractive  index 

spectrum  to  zero  for  those  wavenumbers  which  correspond  to  dimensions 

18 

less  than  the  inner  scale.  Evidence  suggests  that  the  inner  scale  may  be 
typically  between  2. 1  and  5.8  cm  under  neutral  conditions. 

It  is  apparent  from  the  above  considerations  that  the 
extent  of  the  subrange  itself  must  be  carefully  studied.  The  breakdown  of 
the  outer  scale  may  be  observed  with  variable  probe  separations.  The  inner 
scale  is  more  difficult,  since  it  is  not  possible  to  obtain  reliable  measure¬ 
ments  with  very -closely-spaced  probes.  The  use  of  scintillation  experiments 

4 

per  sc  to  measure  the  inner  scale  has  been  suggested,  but  possibly  rep¬ 
resents  a  circular  corroboration  in  the  present  state  of  the  relevant  theories. 
The  most  fruitful  approach  is  thought  to  be  the  measurement  of  the  temporal 
spectrum  of  temperature  fluctuations.^ 
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In  view  of  the  above  discussion,  the  importance  of 

monitoring  the  subrange  --  as  opposed  to  simply  determining  C2  --  is  very 

n 

clear. 


C.  DIMENSIONAL  ASPECTS 

As  discussed  in  the  preceding  report,  the  existence  of  a  break¬ 
point  for  saturation  or  multiple  scattering  raises  serious  dimensional  dif¬ 
ficulties  ,  especially  if  the  point-source  saturated  variance  level  is  definitely 
proven  to  be  wavelength-dependent.  That  is,  if  the  critical  condition  for 
saturation  is  not  simply  that  Eq.  (1)  becomes  some  fraction  of  unity  regard¬ 
less  of  wavelength,  pathlength,  or  C2,  then  the  parameters  of  Eq.(l)  are 

n 

seemingly  insufficient  to  explain  a  breakpoint  from  direct  dimensional  reason- 

25  2 

ing.  Livingston  has  suggested  that  the  single  parameter  is  insuf¬ 
ficient  to  characterize  the  medium  under  the  onset  of  multiple  scattering,  and 
it  is  not  clear  how  the  simple  incorporation  of  a  turbulence  scale  will  alleviate 
the  problem.  This  question  of  a  dimensionless  breakpoint  is  the  single  most 
paramount  issue  in  the  saturation  of  scintillations,  and  in  principle  rises 
above  questions  of  the  validity  of  particular  theories;  all  that  is  needed  is  a 
proper  definition  of  the  relevant  parameters  . 

It  will  be  highly  useful  if  a  Kolmogorov  atmosphere 
can  be  utilized  in  a  substantiated  derivation  of  an  asymptotic  expression  such 
as  that  of  Eq.  (2).  The  result  times  four  may  then  be  set  equal  to  the  right 
side  of  Eq.  (1)  in  order  to  define  a  breakpoint.  For  instance,  a  critical 
distance  L  may  be  defined  as  a  function  of  other  parameters,  and  the  hypo- 

Cr  Vz 

thesis  that  the  amplitude  correlation  length  "saturates"  as  (\  L^)  may  be 
tested . 


The  correct  statement  of  the  dimensionless  break¬ 
points  may  be  somewhat  more  complicated  than  many  dimensional  analyses. 
For  instance,  a  function  which  would  ostensibly  satisfy  the  requirements 
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might  look  like  Eq.  (3): 


(3) 


a 


2 


1 

A 


In  [1  + 


A  oj  i 
(1  +  B  of)  J 


where  o'2  is  the  actual  variance  and  o-2  is  the  Rytov  prediction.  Such  a  func¬ 
tion  gives  <r2  «  o-2  at  low  variances,  exhibits  a  maximum,  and  has  an  asymptote 
(1/A)  In  [l  +  A/B]  .  The  breakpoint  would  be  given  by  cr2  =  (A-B)  1 ,  with  the 
wavelength  dependence,  etc.  appearing  in  the  dimensionless  constants  (A,B). 


IV.  FUTURE  PLANS 


In  view  of  the  above  considerations,  the  program  plan  given  in  the 
preceding  report  may  be  restated  as  follows.  In  order  of  priority,  we  will: 

1.  Gather  comprehensive  point-source  multiwavelength  data 
with  optical  and  thermal  backup  under  saturated 

and  non-saturated  conditions. 

2.  Investigate  the  extent  of  the  inertial  subrange  under 
various  conditions,  as  described  above. 

3.  Investigate  transmitter-aperture  effects  under  saturated 
and  non-saturated  conditions. 

4.  Conduct  receiver  aperture-averaging  experiments. 

The  experimental  systems  and  techniques  for  all  of  these  determina¬ 
tions  are  complete.  The  rate  at  which  the  data  can  be  gathered  will  depend 
primarily  on  weather  conditions. 
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In  order  to  properly  utilize  the  investment  which  has  been  made  on 
this  program,  it  is  suggested  that  a  final,  additional  one-year  period  be 
allotted.  This  will  permit  data -taking  during  the  critical  summer  months 
of  high  turbulence,  and  will  ensure  that  all  of  the  relevant  questions  are  pursued. 
If  time  permits,  this  follow-on  period  could  include  measurements  over  a 
longer  path.  No  further  instrumentation  will  be  developed;  the  entire  effort 
will  be  on  measurements  and  interpretation. 
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APPENDIX  A:  Thermal  Index-of-Refraction  Instrumentation 


The  most  sensitive  means  of  measuring  the  fluctuating  temperature 
difference  between  two  fine -wire  probes  is  to  utilize  an  audio  frequency  car¬ 
rier  driving  a  balanced  bridge  which  contains  the  probes  as  two  of  its  elements. 
This  technique  overcomes  1/f-noise  sensitivity  limitations  and  permits  the 
use  of  a  very-low-noise  front-end  amplifier.  The  tendency  for  the  bridge 
to  continually  drift  out  of  balance  suggests  the  use  of  a  long -time -constant 
feedback-stabilization  feature.  This  system  will  now  be  described  in  detail, 
including  calibration  and  operating  characteristics. 

1 .  Circuit 

The  circuit  is  shown  in  Figure  6.  Modular  (Burr -Brown 
Corp.)  elements  are  utilized  throughout. 

The  bridge  is  driven  with  a  10  kHz  oscillator,  and  its  output 
is  transformer  impedance-matched  to  the  input  of  a  low-noise  amplifier 
(noise  figure  =  2  dB).  After  further  amplification,  and  filtering  in  a  2  kHz 

passive  bandpass  filter,  the  signal  is  fed  to  an  amplifier  which  is  manually 

»,2 

gain-switched  in  six  steps  of  (10)  each.  Following  this,  the  signal  is 
demodulated  in  an  active  precision  rectifier  circuit  and  low  pass  filter.  At 
this  point,  an  output  is  provided  which  is  linearly  proportional  to  the  tem¬ 
perature  fluctuations  (AT),  for  purposes  of  probability  and  spectral  analysis. 

The  demodulated  signal  (AT)  is  then  a. c. -coupled  with  a 
long  (0.8  sec)  time  constant  through  buffers  to  a  bipolar  squaring  module 
which  current-feeds  the  summing  point  of  a  following  operational  amplifier. 

In  the  final  stage,  the  resulting  d.c.  signal  <AT*>  is  averaged  in  an  integrator 
which  is  provided  with  short  (1.8  sec)  and  long  (1  minute)  time  constants. 

The  result  appears  on  a  500  microamp  meter  and  is  available  as  a  d.c.  out¬ 
put  voltage. 
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The  bridge  constitutes  a  "balanced  modulator"  with  cancelled 
carrier,  which  would  result  in  an  output  envelope  proportional  to  |  AT  |  .  To 
restore  the  carrier  and  hence  achieve  the  desired  linear  amplitude  modula¬ 
tion,  it  is  possible  to  use  gated  or  synchronous  detection  with  a  reference 
input  from  the  oscillator.  However,  a  small,  deliberate  bridge  imbalance 
creates  the  same  effect  while  providing  a  simple  means  of  achieving  feed¬ 
back  stabilization. 

The  output  of  the  demodulator -filter  is  d.c.  coupled  through 
a  slow  integrator  to  optically- coupled  variable-resistor  elements  which  are 
in  parallel  with  the  bridge  arms.  The  degree  of  bridge  imbalance  is  con¬ 
trolled  by  an  adjustable  d.c.  offset  at  the  input  to  the  integrator.  The  loop 
gain  is  maintained  constant  as  the  signal  gain  is  varied,  by  means  of  a  coupled 
switch  as  shown. 

The  action  of  the  feedback  circuit  is  much  slower  than  the  a.c. 
coupling  time  constant  in  the  output  signal  path,  so  that  slow  drifts  are  can¬ 
celled  without  affecting  the  reading.  This  is  permissible  because  the  static 
temperature  or  resistance  difference  between  the  probes  is  not  of  interest. 

In  order  to  provide  a  small,  incremental  bridge  imbalance 
free  of  hand  effects,  for  calibrating  purposes,  a  reed  relay  is  utilized. 

When  driven  by  a  slow  square-wave  generator,  this  feature  permits  easy 
calibration  and  a  gain- stability  check  of  the  system. 

In  order  to  realize  the  maximum  possible  dynamic  range,  the 
levels  are  chosen  so  that  saturation  (±  10  volts)  occurs  simultaneously  in 
various  elements.  Instantaneous  saturation  peaks  are  indicated  by  means 
of  a  peak  detection  circuit  and  light- emitting  diode. 

The  circuit  is  powered  with  rechargeable  Ni:Cd  cells.  Meter 
switching  at  the  output  provides  capability  for  monitoring  the  battery  levels, 
carrier  (bridge  offset)  level,  and  output  proportional  to  C  2  . 
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2. 


Operation 


Operation  of  the  system  is  quite  simple.  With  the  feedback 
loop  opened  up  and  caps  on  the  probes,  the  bridge  is  balanced  by  monitoring 
the  demodulator  output  and  manually  trimming  the  resistance  and  capacitance 
at  the  bridge.  The  feedback  is  engaged  and  set  for  a  carrier  level  equal  to 
half  the  saturation  voltage,  to  permit  a  maximum  linear  modulation  range. 
(This  carrier  condition  is  maintained  automatically  as  the  gain  is  varied.) 

The  unit  is  faced  into  the  wind  with  the  desired  height  and  probe  separation, 
and  the  probes  are  uncapped.  The  gain  is  increased  until  occasional  satura¬ 
tion  is  noted,  and  then  backed-off  one  step.  Each  gain  step,  after  squaring, 
represents  a  change  in  sensitivity  by  a  factor  of  ten. 

The  system  will  stably  monitor  at  the  output  meter.  It  is 
necessary  to  record  the  meter  reading,  gain  setting,  barometric  pressure, 
and  temperature.  Through  the  use  of  a  long  coaxial  cable  and  the  scintilla¬ 
tion  computer  electronics,  the  probability  distribution  and  spectrum  of  the 
temperature  fluctuations  are  readily  analyzed;  in  the  case  of  the  spectrum, 
one  probe  may  be  capped  so  that  absolute  fluctuations  are  measured. 

If  the  unit  is  taken  from  a  warm  environment  into  a  cold 
weather  condition,  the  capacitance  drift  may  eventually  require  rebalancing. 
In  any  fully-automatic  version  of  this  system,  such  as  in  balloon  experi¬ 
ments,  it  will  be  necessary  to  use  temperature-compensated  circuitry,  or 
alternatively,  quadrature  synchronous  detection  may  be  employed  with 
varactor  diodes  to  automatically  cancel  capacitance  imbalance. 

3 .  Linearity  Considerations 

The  bridge  output  is  related  to  an  incremental  imbalance  by 

a  simple  Thevenin  equivalent  circuit.  If  we  let  Rj  =  fixed  bridge  arms,  R2 

and  (R2  +  AR)  =  the  probe  resistances,  R  =  oscillator  source  resistance, 

s 
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V  -  oscillator  voltage,  and  R  -  transformer -coupled  load  resistance,  the 
8  L 

bridge  output  may  be  shown  to  be 


V  =  V 
out 


s 


RtAR  Rl 


<Ri  + 


R2)  (2Rs+ 


Ri+  R2) 


(R  +  2  RtR2) 
Ri+R2 


(A- 1) 


In  this  expression,  it  is  assumed  that  AR/R2«  1.  Now,  AR  will  consist 
of  a  deliberate  offset  or  "bias",  plus  the  a.c.  temperature  fluctuations. 

Under  conditions  of  highest  turbulence,  the  thermal/resistance  coefficient 
of  the  platinum  probe  wires  is  such  that  AR/R2  is  on  the  order  of  10  2,  which 
also  determines  the  maximum  necessary  bias.  Hence,  linearity  is  not  notice 
ably  degraded. 


The  effects  of  capacitive  imbalance  may  be  examined  by 
replacing  AR  above  by  AZ .  (Symmetrical  or  balanced  capacitance  components 
will  also  affect  the  bridge  gain;  e.g.  ,  replace  Ri  by  Zj).  It  is  found  that 
linearity  is  maintained  if  cjR2AC  «  AR^  /R2,  where  w  is  the  carrier  radian 
frequency.  This  is  only  a  problem  for  long  unattended  periods  under  lowest 
turbulence  conditions,  where  AR^.^  is  small.  Automatic  capacitance  trim 
is  readily  employed,  as  discussed  above,  but  is  not  necessary  under  our 
experimental  conditions. 


4.  Sensitivity  Considerations 

The  system  described  here  nearly  achieves  the  fundamental 
sensitivity  limits  set  by  Johnson  noise.  The  fundamental  sensitivity  of  the 
system  may  be  calculated  by  first  assuming  a  noiseless  front-end  amplifier, 
which  has  an  equivalent  input  rms  noise  voltage  given  by 

1  Vz 

V  =  7  (4KTBR.)  .  (A-2) 

n  c  l 
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In  this  expression,  T  may  be  taken  as  300°K,  B  is  2  kHz  (carrier  band¬ 
width),  and  is  the  amplifier  input  resistance.  This  may  be  related  to 
an  equivalent  AT  at  the  probes.  For  simplicity  (and  to  realize  the  best  noise 
figure),  we  take  =  R.  The  oscillator  current  is  Ifl  =  V^/R^,  with 

R  »  R  so  that  a  current  source  may  be  assumed.  If  the  transformer  turns 
s  i 

ratio  is  N  =  (R^/RjJ  »  we  have 

I 

V  .  =  ~  AR  N  .  (A-3) 

sig  rms  8  rms 


The  SNR  is  given  by 


I  2 

SNR  =  r-2-  (AR)2  N2/kTB  N^ 

64 


I  2  R 

8 

64kTB 


(A -4) 


P  a2  (AT)2 
16  kTB 


(A -5) 


where  a  is  the  temperature  coefficient  of  the  probe  wire  (AR/AT),  and  P 

is  the  probe  heating  power  (I  2/4R).  The  minimum  or  critical  temperature 

s 

sensitivity  is  thus 


AT 


cr 


4  (kTB) 


72 


a  P 


/2 


(A-6) 


The  value  of  P  is  limited  by  permissible  probe  heating,  since  sensitivity 
to  wind  speed  is  undesirable.  For  the  present  case,  I  =  2.4  x  10~4,  P  =  10~6, 
a  =  3  x  10  3,  and  hence  ATcr  =  3.2  x  10  3  “C.  This  figure  must  be  in¬ 
creased  according  to  the  noise  figure  of  the  amplifier,  which  in  our  case  is 
2  dB.  This  is  confirmed  by  direct  measurment. 
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The  above  appraisal  is  overly  conservative  when  it  is  realized 
that  noise,  within  the  2  kHz  bandwidth,  merely  contributes  an  offset  in  the 
squared,  averaged  output.  The  constancy  of  this  offset  is  what  really  counts, 
and  hence  the  effective  noise  bandwidth  is  on  the  order  of  the  reciprocal  of 
the  integration  time,  or  4-5  orders  of  magnitude  less  than  the  2  kHz  value 
used  above.  For  example,  with  the  present  unit  set  at  its  highest  gain,  with 
which  one  may  readily  measure  turbulence  levels  lower  than  ever  realized 
at  the  ground  (see  below),  the  noise-offset  is  still  not  discernible  on  the  meter. 

The  choice  of  probe  wire  diameter  is  an  important  design 
factor.  This  choice  is  a  trade-off  between  fast  response  on  one  hand,  and 
ruggedness  and  higher  sensitivity  on  the  other.  The  higher  sensitivity  with 
larger  wires  is  due  to  the  larger  probe  heating  power  P  which  is  permitted 
as  the  surface  area  increases.  The  response  to  fast  AT  components  can  be 
increased  electronically,  but  noise  is  also  enhanced.  In  the  present  system, 
chopped-laser  measurements  showed  adequate  thermal  response  with  probes 
of  2-micron  diameter,  where  P  is  on  the  order  of  1  microwatt. 


5.  Calibration  and  Measurement  Range 

The  meter  output  is  related  to  the  square  of  the  fractional 
bridge  imbalance,  as  follows: 


( AR) 

V  =  K  -4~- 
out  R2 


(A -7) 


which  defines  K.  The  temperature  structure  constant  ^  is 


C  2  = 
T 


<AT2>  (AR/qR)2 
r2/3  = 


out 

r2'3  K  a  z 


(A-8) 
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where  r  is  the  probe  separation.  In  this  expression,  the  factor  K  is  deter¬ 
mined  through  the  use  of  the  calibrating  AR  built  into  the  instrument,  as 
described  above.  The  final  necessary  expression  relates  C  2  to  CT 
(Refs.  12-14): 

C  =  ■7.7-A-P.  (i  +  0. 00753 /V2)  10~6  C_,  (A-9) 

n  T‘  i 

where  P  is  barometric  pressure  in  millibars,  X.  is  in  microns,  and  T  is  in 
°K.  For  operating  convenience,  this  is  reduced  to  an  expression  at  each 
wavelength  in  terms  of  the  pressure  in  inches  and  the  temperature  in  °F. 

The  parameter  K  is  a  function  of  the  gain  setting.  The  useable 
range  of  the  instrument  covers  values  of  from  1.24  x  10  11  m  Z^J  to 
1 . 24  x  10~17m"z/3,  which  more  than  brackets  values  ever  observed  at  the 
ground. 


6 .  Spectral  and  Probability  Analysis 

As  with  the  scintillations,  the  spectrum  of  the  AT  signal  is 
analyzed  with  a  Tektronix  Type  1L5  unit.  Since  this  instrument  is  linear  in 
amplitude,  a  falloff  exponent  of  -5/6  vs.  frequency  is  expected  according  to 
the  one -dimensional  Kolmogorov  result.  To  examine  this  in  detail,  it  is 
possible  to  use  the  analyzer  in  a  logarithmic  mode,  so  that  the  slope  of  a 
log-log  plot  may  be  precisely  determined. 

The  probability  distribution  of  AT  is  analyzed  in  the  scintil¬ 
lation  computer.  AT  is  related  to  the  linear  output  voltage  through  gain 
calibration,  and  the  variance  <ATZ>  is  determined  from  the  slope  on  gaussian 
probability  paper.  is  then  determined  from  (A-8)  and  (A-9).  The  presence 

of  nonlinear  "tails"  on  the  probability  plot  indicates  components  outside  the 
inertial  subrange,  as  discussed  in  Section  III. 
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APPENDIX  B:  Notation,  Definitions,  and  Relationships 
for  Scintillation  Quantities 

1 .  Basic  Quantities 

u  =  complex  optical  wave  amplitude  (rms) 

u  =  <  lul  > 

o 

l  =  In  ( |uj  /u^)  =  log  amplitude  (normalized) 

4>i  =  In  (u/u  )  =  Rytov  wave  function 
o 

t  =  Re  (4*1) 

I  =  intensity  =  |u| 2 

I  =  u  2 
o  o 

f  =  1/2  In  I/I 

o 

Important  identity:  i  -  <i>  =  ln|u|  -  In  u  -  <lnlul>  +  <ln  u  > 

o  o 

=  ln(ul  -  <lnlul> 

Hence,  the  log  amplitude  is  obtained  electronically 
by  a.c.  -  coupling  the  logarithmic  signal. 


2 .  Moments: 

C^(0)  =  <(i  -<f>)2>  -  log  amplitude  variance 
Cj(0)  =  <(ln  I  -  <ln  I>)2  >  =log  intensity  variance 
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(I  -  <I>)2 

<I>2 

<I2> 


In 


<I>2 


normalized  intensity  variance 


3 .  Relationships; 


CjfO)  =  4^(0) 
azT  =  In  (1  +  <r2  ) 


If  the  statistics  are  log  normal,  then 


'i  =  e 


CjIO) 


Cj(0)  =  cr2T 
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R  AC  T 

In  conducting  multiwavelength  scintillation  experiments  over  a  one 
mile  path,  it  was  found  that  a  finite-sized  transmitter  aperture  does 
not  represent  a  true  spherical-wave  source,  regardless  of  the  beam 
divergence  used.  This  observation  has  important  implications  for  the 
interpretation  of  other  data  on  the  "saturation"  of  scintillations, 
and  in  the  present  experiments  has  lead  to  the  redesign  of  the  transmitter 
to  achieve  a  virtual  point  source.  This  has  been  completed,  including  a 
provision  for  independently  increasing  the  aperture  at  each  wavelength 
up  to  two  inches  in  diameter,  for  later  transmitter-aperture-dependence 
measurements . 

A  highly-developed  thermal  probe  system  for  the  detailed  6tudy  of 
turbulence  structure  was  made  completely  operational,  and  a  portable 
laser  system  was  fabricated  for  independent  "optical  turbulence  strength" 
measurements.  All  subsystems  of  the  experimental  facility  are  operating 
without  problems,  and  no  further  modifications  are  contemplated. 

In  order  to  put  the  problem  of  atmospheric  scintillations  in  current 
perspective  and  to  guide  the  comprehensive  experiments,  a  discussion  of 
recent  experimental  and  tneoretical  issues  is  given.  It  is  felt  that 
the  understanding  of  the  turbulence  scattering  problem  has  progressed 
to  the  point  where  the  correct  questions  may  be  posed  and  definitively 
investigated  through  these  experiments. 
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